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Abstract 



We calculate nuclear enhanced power corrections to structure functions mea- 
sured in deeply inelastic lepton-nucleus scattering in Quantum Chromody- 
namics (QCD). We find that the nuclear medium enhanced power corrections 
at order of 0(a s /Q 2 ) enhance the longitudinal structure function F^, and 
suppress the transverse structure function F\. We demonstrate that strong 
nuclear effects in oaI^d and Ra/Rd-, recently observed by HERMES Collab- 
oration, can be explained in terms of the nuclear enhanced power corrections. 
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Recently, a strong nuclear dependence in lepton-nucleus deeply inelastic scattering (DIS) 
was observed by HERMES collaboration Q. HERMES data show a large nuclear enhance- 
ment in the longitudinal cross section when Bjorken xb is small. The ratio of longitudinal 
cross sections cr^/crf measured by HERMES Collaboration was as large as 2 for scattering 
off a 14 N target. On the other hand, due to nuclear shadowing in the small x region, the 
ratio of the transverse cross section a^/a^ is much less than 1, which leads to a very large 
nuclear enhancement in a double ratio / a T) / ( a L /°r) ED- The anomalous nuclear ef- 
fects observed by HERMES have generated lots of interests p|. In this Letter, we calculate 
the nuclear dependence of the DIS cross sections (or DIS structure functions) at low Q 2 in 
terms of medium enhanced power corrections in QCD perturbation theory, and provide our 
understanding of HERMES data. 

For low energy fixed target DIS processes, like HERMES experiment, we neglect a small 
contribution from vector boson Z and take the one-photon exchange approximation. Like 
other DIS experiments, HERMES experiment measures inclusive lepton-nucleus cross sec- 
tions as a function of xb and virtual photon's invariant mass square Q 2 JH 

dojh q~Mott 7iF 2 (x B ,Q 2 ) 
dxsdQ 2 EE' e xb 

where OMott represents an elastic cross section for lepton scattering from a point charge, and 
E and E' are incoming and scattered lepton energy, respectively. The parameter e is given 
by! 

A{l-y)-Q 2 /E 2 
e 4(1 - y) + 2y 2 + Q 2 /E 2 ' 

where y = (E — E')/E. The R(xb,Q 2 ) in Eq. ([]]) is the ratio of longitudinal to transverse 
DIS cross sections R = a^jaT- When xb < 0.03, HERMES data on the ratio of inclusive 
DIS cross sections from nucleus A and deuterium D(= 2 R), cta/cd, show major differences 
from what have been measured by other experiments in terms of the dependences on xb, 
Q 2 , and e 0. When xb decreases, the ratio falls steeply, and is much smaller than what 
was measured by other experiments. When Q 2 increases, the ratio is not only small but also 
decreases. On the other hand, the ratio increases when e increases. Why HERMES data is 
so different from other data in the small xb region? 

From Eq. (|I]), the ratio of inclusive DIS cross sections from nucleus A and deuterium 
D(= 2 H) can be expressed as 

a A _ F 2 A (l + eR A )(l+R D ) 
a D F 2 D (l + R A )(l + € R D ) ' 

where Ra and Rd represent the ratio of ctl/ctt for nucleus A and deuterium D. HERMES 
Collaboration shows that the ratio of F 2 / F 2 D extracted from its data is consistent with the 
previous measurements, and attributes the differences in the ratio of cross sections, (Ja/^d, 
to the strong nuclear dependence in Ra/Rd [0]- At the same xb, the key difference between 
HERMES data and other data is the range of Q 2 and e. When xb is small, the Q 2 of 
HERMES data is much smaller than that of other data sets. It is then natural to investigate 
the role of 1/Q 2 type power corrections to the ratio of the DIS cross sections, R. 
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l + eR{x B ,Q 2 ) 
l + R(x B ,Q 2 ) 



(2) 



(3) 



When energy exchange of the collision, Q, is much larger than the typical momentum 
scale of the hadron wave function, the DIS process is dominated by a single hard collision 
between the lepton and a quark via the virtual photon, as shown in Fig. |l|(a). The cor- 
responding cross section is given by the probability to find a quark within a hadron times 
an elastic cross section between this quark and the incoming lepton. QCD radioactive cor- 
rections to such a single parton scattering picture are well-defined 0. However, in a large 
nucleus, partons from different nucleons can interact with each other before the hard colli- 
sion takes place, as shown in Fig. 0(b), or interact with the scattered quark after the hard 
collision, as shown in Fig. [TJ(c). Just like the process in Fig. [TJ(a), leading contributions from 
the process shown in Fig. |l](b) is dominated by the phase space where the quark of momen- 
tum k is pinched to be on-mass shell. That is, the multiparton interactions from different 
nucleons are separated from the hard collision by a "long" time scale. Therefore, such mul- 
tiparton interactions are internal to the nucleus, and do not affect the short- distance hard 
collision between the quark and the virtual photon. On the other hand, such interactions 
are responsible for the differences between the parton distribution in a nucleus and that in 
a free nucleon |4]] . Because of the leading twist nature of the operators defining the parton 
distributions in a nucleus, and the fact that F 2 is dominated by the contributions from the 
parton distributions, the ratio of F 2 / F 2 has a small Q 2 dependence [§]•§[. 

Unlike the parton momentum k in Fig. |l|(b), the parton momentum k' in Fig. [l|(c) is 
not pinched after we square the scattering amplitude due to different gluon momenta in 
the amplitude and its complex conjugate 0. The integration of the unpinched parton 
momentum k' can be deformed away from the k' 2 = 0, and the quark-gluon interaction in 
Fig. |l](c) is a part of the hard collision. Because of the extra interactions with the scattered 
quark, physical contributions from such multiparton scattering are suppressed by inverse 
powers of Q 2 , known as the high twist contributions P-TT0|]. Although their contributions 
to the DIS structure functions are normally small in lepton-hadron collisions, such power 
corrections can be important in the nuclear collisions because of the enhancement of nuclear 

H 



size 



In the rest of this Letter, we show that the power corrections caused by the type 
of multiparton interactions shown in Fig. |l|(c) are responsible for the features seen in the 
HERMES data. 

In terms of DIS structure functions, F 1 (xb,Q 2 ) and F 2 (xb,Q 2 ), the ratios of the DIS 
cross sections can be expressed as |12j 



<Jl_ 



2x B F l 



1 + 



AmijX 



Q 2 



F 2 - 2x B F 1 



2F, 



(4) 



where is nucleon mass, and the longitudinal structure function Fl is defined as [] 



F L = — 

x b 



^m 2 N x 2 B 



Q 2 



F 2 - 2x B F 1 



— [F 2 - 2x B F 1 ] 
x B 



(5) 



when xb "C 1 or m 2 N <C Q 2 . Let us introduce two projection operators, 



x Our definition of Fl in Eq. (||) differs by a factor of 1/xb from that in Ref. [jij. 
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cf = ^ W + ZxbjT] [q v + 2x B p v \ , (6) 

1 2r 

eg? = [p/y + 9 V] + — - fl r ' 11/ , (7) 

p ■ q p ■ q 

where is virtual photon's four-momentum with q 2 = —Q 2 . The DIS hadronic tensor W^ v 
can be then expressed as || 

W^(x B ,Q 2 ) = ~ [ef F L (x B ,Q 2 ) +e^2F 1 (x B ,Q 2 )] . (8) 

Therefore, we can extract the DIS structure functions from the hadronic tensor as 

F L (x B ,Q 2 ) = 2e1 u W, u (x B ,Q 2 ), and F x (x B , Q 2 ) = ^ e^W iiV [x B , Q 2 ) . (9) 

At the lowest order, the DIS structure functions in lepton-nucleus collisions are given by 

F£(x B ,Q 2 ) = 0, (10) 
Ff{x B , Q 2 ) = ^ F 2 A (x B , Q 2 ) = \Y. 4 <f> g /A(xB, Q 2 ) , (11) 

where J2 q runs over all quark and antiquark flavors and e q is corresponding fractional charge. 
The (p q /A is an effective quark distribution in a nucleus of atomic number A, which is nor- 
malized by dividing the A. 

The lowest order power corrections to the DIS structure functions are given by the 
Feynman diagrams in Fig. |2|, plus the diagrams with four-quark lines instead of two-quark- 
two-gluon lines |§. Since we are interested in the medium size enhanced power corrections, 
we neglect the four-quark processes [11]. In Fig. |^, the quark lines with a short bar are 
special quark propagators, which are equal to the normal quark propagators with the leading 
twist contributions removed 0. The diagrams with special quark propagators in Fig. |2| 
are important to preserve the electromagnetic current conservation at the level of power 
corrections. Feynman rules for the special quark propagators are given in Refs. P,|T3||. In 
the light-cone gauge of the strong interaction, contributions to the hadronic tensor from the 
diagrams in Fig. ^] can be factorized as [ IT ] 

WjjX(x B ,Q 2 ) = — J dxdxidx 2 T qA (x,Xi,x 2 )H^(x, x 1 ,x 2 ,q) • (12) 
The hadronic matrix element T qA is given by 

T qA (x, Xl ,x 2 ) = ( ( f^l^L^l e ^P + y~ e ^) P + yT e -H^)P^ 
\x — X\J \x — x 2 J J 2n 2n 2n 

x (P A \m^F +a {y 2 )F+{y^m V -)\P A ) , (13) 

where p = P A /A is averaged momentum per nucleon. In deriving T qA , we used F +a (y~) = 
n p d p A a (y~) in the light-cone gauge. The partonic part H^ A is given by the Feynman dia- 
grams in Fig. § with the quark lines contracted by ^7 -p and gluon lines contracted by \dap- 
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The transverse tensor d a p = —g a p + n a np + n a np with two lightlike vectors, n M = and 
n ti = After contracting with the projection operators, we obtain || 

e£H#(x B , Q 2 ) = 4 e\ { i *(* - x B ) , (14) 



e^<(x B ,Q 2 ) = 2e;(^j-x B 




<5(x 2 — xb) S(xi — xb) 



X\ — Xb X2 — Xb 



(15) 



By substituting the above partonic hard parts into Eq. flI2|), we can derive the leading power 
corrections to the DIS structure functions defined in Eq. @. 

In order to derive the medium size enhanced power corrections, we need to identify the 
position space integrals that can be extended freely to the size of the target. For the leading 
power corrections in Eq. fll2|) , only position space integrals are from the dy~'s in Eq.(|i~3|). 
Generally, these y~ integrals cannot grow with the size of the target because of oscillation 
of the exponential factors. Because of the nature of Fourier transform, the exponentials 
in Eq. (|13| ) are linear functions of parton momentum fractions. Since the kinematics of 
inclusive DIS is only sensitive to the total momentum from the target, two of the three 
momentum fractions in Eq. (|T2|) cannot be fixed by the hard collision. As pointed out in 



Ref . Jl 1| , the integration of these two momentum fractions can be fixed by two unpinched 
poles from the partonic parts, and the integration removes the oscillation of two exponential 
factors. Therefore, two corresponding y^ integrations can be extended freely to the whole 
target, and generate the target size enhanced power corrections. 

After combining the hadronic matrix element in Eq. ([13]) with the partonic hard parts in 
Eq. flT4]), the momentum fraction x is fixed by the S(x — ig), while dx\ and dx-i integrations 
are dominated by the region where X\ ~ xb and x 2 ~ xb, because of the two poles at 
l/O^s — ^i) an d ]-/{xb ~ X2) respectively. The fact that two poles are from the matrix 
element and do not have explicit ie is an artifact of choosing the light-cone gauge. In order 
to recover the correct ie within the light-cone gauge calculations, we adapt the ie prescription 
introduced in Ref. [Q, 

* ' 1 W 1 ^ (16) 



,Xb — XxJ \Xb—X2/ \Xi—Xb+1£/ \X2 — xb — it, 

We can also work in a covariant gauge to avoid the spurious poles in the light-cone gauge 



T3|,n]. After carrying out dx\ and dx 2 integrations, we obtain the leading power corrections 
to the longitudinal structure function, 

F l^Q\ /q2 = £ej4 ^T> B ,Q 2 ) , (17) 

where J2 q runs over all quark and antiquark flavors, and the twist-4 quark-gluon correlation 
function is given by [|lj|] 

T q A F (x B , Q 2 ) = J d -t elXBP+y ~ J dj ^r- %r - v-) % 2 ~) 

x (pa\M0) X-F +a ( y2 -)F a + (y^ q (y-)\PA) ■ (18) 
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In order to derive leading power corrections to Fi, we expand the x\ and £2 in Eq. (p75|) 
around x because of the poles in Eq. (|13"D , 



5(x 2 — x B ) S(xi — x b 



+ 



-8'(x — x B ) + 



(19) 



x\ — xb X2 — xb 

where represents terms proportional to powers of Xj — x with % = 1,2. Since these 

terms cancel the poles in Eq. (|T3|), we can neglect them for calculating the leading nuclear 
size enhanced contributions, and we obtain 



Ff(x B ,Q 2 



i/Q 2 




d 



■Ti 



dxB -A x B,Q 



(20) 



where J2 q runs over & U quark and antiquark flavors, and T A F is given in Eq. (|i~8|). 

In order to estimate the size and sign of the leading power corrections to the structure 
functions in Eqs. ( |TTD and (PO]), we need the knowledge of the twist-4 quark-gluon correlation 
function T a f {xb, Q 2 )- Although this function is nonperturbative, QCD factorization theorem 
enables us to find the information on T A F from other physical observables. For the following 



numerical estimate, we adopt the model [15 



F^f(x, Q 2 ) 



X 2 A l ^<f> q/A {x,Q' 



(21) 



and use A 2 

in hadron-nucleus collisions 



0.01 GeV 2 extracted from the Drell-Yan transverse momentum broadening 
Substituting T A F in Eq. 



TBI 



into Eqs. ([T7|) and (gDJ), we 
find that the leading nuclear size enhanced power corrections to Fl are positive, while the 
corrections to F\ are negative. Because of the opposite signs, the nuclear enhanced power 
corrections to F2 are smaller. 



In order to understand the HERMES data, we consider the following double ratio, 

[l + eR A ) /(l + eR D ) 



Rr = 



[1+Ra 




l + R 



D 



(22) 



Since F A /F 2 D measured by HERMES Collaboration is consistent with the ratio measured 
by other experiments Q, the double ratio Rr should be responsible for all unique features 
of HERMES data, such as the xb, Q 2 , and e dependence of the ratio oaI^d-, defined in 
Eq. (0). From Eq. (|), we have R A = F£/(2Ff) with the power corrections included for the 
DIS structure functions 



F t A (x B ,Q 2 ) = Ff(x B ,Q 2 ) T +F l A (x B ,Q 



LT 



1/Q 2 



(23) 



where i — L, 1 and "LT" represents the leading twist contributions. The leading power 
corrections to Fi and F\ are given in Eqs. ( |TTD and (pO"D, respectively. The leading order 
leading twist contributions are given in Eqs. (0) and (|10D, respectively. Since the leading 
order F A = and the order of a s contributions are renormalization scheme independent, we 
use them as the nonvanishing contributions to Fl (T7 



[I 


mi 











l xb 
x' 



(24) 
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Since Rd is well measured and Rd ~ Rp of proton |fL8[j , we let Rd = F£ /{2F^) with 
only the leading twist contributions. We find that our Rd in low Q 2 and low xb region are 



consistent with the parameterization in Ref. ||18||. Because we do not have well-tested parton 



distributions at Q 2 as low as 0.5 GeV 2 , we will use a set of generic parton distributions in 
Ref. [|19j for the following numerical estimate. 

With the model of T^ F in Eq. (|2l|) and A 2 = 0.01 GeV 2 , and without any further free 
parameter, we plot the xb and Q 2 dependence of the double ratio Rr for 14 N in Fig. 0(a) 
and 0(b), respectively. In plotting Fig. 0, we used E = 27.5 GeV for the lepton beam 
energy. In Fig. ^(a), the ratio is sensitive to the value of Q 2 , and Q 2 = 0.6 GeV 2 was 
used. In Fig. 0(b), three curves correspond to xb = 0.0175 (dotted), 0.025 (solid), and 0.035 
(dashed), respectively. Clearly, Figs. 0(a) and 0(b) reproduce all the key features of the 
HERMES data shown in Fig. 1 and 2 of Ref. (I| . The detailed numerical comparison with 
the HERMES data requires better knowledge of original data analysis and reliable low Q 2 
parton distributions, and will be presented elsewhere. 

In the rest of the Letter, we explain why the nuclear enhanced power corrections are 
consistent with the HERMES data. As shown in Fig. 0(a) or Fig. 1 of Ref. the ratio Rr 
or cta/vd shows a steep falling feature when xb decreases, which is not observed in other 
experiments. We find that the steep falling was mainly caused by strong xb dependence 
of the e and a large Ra- Substituting y = Q 2 j{2m^XBE) into Eq. (0), we find that e falls 
steeply from 0.9 to 0.25 when xb decreases from 0.03 to 0.0125 at Q 2 = 0.6 GeV 2 . Thus, 
the ratio (1 + €Ra)/(1 + Ra) falls along with the e. 

The Q 2 dependence, as shown in Fig. 0(b) or in Fig. 2 of Ref. [l[], is the most counterin- 
tuitive. Normally, when Q 2 increases, we expect the effect of power corrections to go away, 
or the ratio Rr to approach one. However, the HERMES data shows an opposite behavior. 
We find that such a counterintuitive behavior in HERMES data was caused by the strong Q 2 
dependence of the e parameter. For example, when Q 2 changes from 0.5 GeV 2 to 0.9 GeV 2 , 
the e at xb = 0.0175 decreases from 0.78 to 0.13 by dropping a factor of 6, while the size of 
power corrections changes for less than a factor of 2. We find that the decrease of the double 
ratio Rr when Q 2 increases is due to the steep falling of the combination e Ra in Eq. (p2|) . 

Similarly, the e dependence of the HERMES data, shown in Fig. 3 of Ref. [J], is a direct 
consequence of explicit e dependence in Eq. (0) or Eq. (0). From our numerical estimate, 
we also find that the observed large double ratio Ra/Rd ~ 5 for 0.01 < x < 0.03, as shown 
in Fig. 5 of Ref. 0, is due to a strongly nuclear enhanced Ra ~ 1-6 at Q 2 = 0.6 GeV 2 and 
a relatively small Rd ~ 0.32, which is consistent with other measurements P"8"[ . 

We conclude that strong nuclear effects in (JaI^d and Ra/Rd, observed by HERMES 
Collaboration, are direct consequences of strong nuclear enhanced power corrections. We 
demonstrated that our numerical estimate reproduces all the features of the HERMES data. 
When the Q 2 is as small as 0.6 GeV 2 , we find that the structure function Fi(xr,Q 2 ) is 
still dominated by the leading twist contributions, while the longitudinal structure function 
Fl{xb,Q 2 ) is dominated by the nuclear enhanced power corrections. Since the leading 
power corrections are large, it is therefore important to study the nuclear enhanced power 
corrections at even higher powers of l/Q 2 to have a better QCD prediction. 
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FIGURES 




FIG. 1. Sketch for lepton-hadron and lepton- nucleus DIS: (a) single-parton scattering, (b) 
multiparton interaction internal to the nucleus, and (c) multiparton interaction at the hard collision. 




FIG. 2. Leading order Feynman diagrams that contribute to the leading power corrections to 
the DIS structure functions. 
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FIG. 3. The double ratio Rr in Eq. (f22^) as a function of xb (a) and Q 2 (b). The curves are 
explained in the text. 
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